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Abstract We discuss recent results and future prospects regarding the investiga- 
tion, by lattice simulations, of the non-perturbative properties of QCD and of its 
phase diagram in presence of magnetic or chromomagnetic background fields. Af- 
ter a brief introduction to the formulation of lattice QCD in presence of external 
fields, we focus on studies regarding the effects of external fields on chiral sym- 
' . metry breaking, on its restoration at finite temperature and on deconfinement. We 

qj* conclude with a few comments regarding the effects of electromagnetic background 

<**] . fields on gluodynamics. 



1 Introduction 

The properties of strong interactions in presence of strong external fields are of great 
phenomenological and theoretical importance. Large magnetic or chromomagnetic 
background fields, of the order of 10 16 Tesla, i.e. yjcjl? ~ 1.5 GeV, may have been 
produced at the time of the cosmological electroweak phase transition [T) and may 
have influenced the subsequent evolution of the Universe, in particular the transition 
from deconfined to confined strongly interacting matter. Slightly lower magnetic 
fields are expected to be produced in non-central heavy ion collisions (up to ~ 10 15 
Tesla at LHC EH), where they may give rise, in presence of non-trivial topological 
vacuum fluctuations, to CP-odd effects consisting in the separation of electric charge 
along the direction of the background magnetic field, the so-called chiral magnetic 
effect EJ |2j 15). Finally, large magnetic (or even chromomagnetic J6|) fields are ex- 
pected in some compact astrophysical objects, such as magnetars Q (see O for a 
recent review on the subject). 
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Apart from the phenomenological issues above, there is great interest also at a 
purely theoretical level, since background fields can be useful probes to ge insight 
into the properties of strong interactions and the non-perturbative structure of the 
QCD vacuum. A typical example is chiral symmetry breaking, which is predicted 
to be enhanced by the presence of a magnetic background, a phenomenon which is 
known as magnetic catalysis [E3 [HI [12 03 DICE [IS El 

|^^^|^l|2Sl^ll2QllIIll2ll22ll24ll2lll23l22l- 

Such interest justifies the large efforts which have been dedicated in the recent 
past to this subject by a variety of different approaches, most of which are reviewed 
in the present volume. Lattice QCD represents, in general, the ideal tool for a first 
principle investigation of the non-perturbative properties of strong interactions. This 
is particularly true for the study of QCD in presence of a strong external magnetic or 
chromomagnetic field, since no kind of technical problem, such as a sign problem, 
appears, preventing standard Monte-Carlo simulations, as it happens instead at finite 
baryon chemical potential. 

Similarly to the continuum theory, the magnetic field is introduced on the lat- 
tice in terms of additional 1/(1) degrees of freedom (see Sec. 12. U . which are not 
directly coupled to the original SI/ (3) link variables and affect quark propagation 
by a modification of the covariant derivative (i.e. the fermion matrix). The fermion 
determinant, contrary to the case of a baryon chemical potential or of an electric 
background field, is still real and positive, allowing for a probabilistic interpretation 
of the path integral measure. Lattice studies including the presence of electromag- 
netic fields have been done since long, originally with the purpose of studying the 
electromagnetic properties of hadrons f38l l39l l40l I4T1 l42l . The introduction of a 
chromomagnetic field requires a different approach, since in this case the back- 
ground variables are strictly related to the quantum gluon degrees of freedom: a 
standard procedure is that defined in the framework of the lattice Schrodinger func- 
tional E3lll4l|35l|46]|l7l|48]|l9l|50l[5ll (see Sec.E2). 

The recent interest about external field effects in the QCD vacuum and around 
the deconfining transition has stimulated lattice investigations by many groups ll52l 
EH |54l [551 ESI |57l [581 HU [601 HI] [62 |63] HU [651 US). Several studies have fol- 
lowed a quenched or partially quenched approach, considering only the effect of 
the magnetic background on physical observables, i.e. on quarks propagating in the 
background of the given magnetic field and of non-Abelian gauge configurations, 
the latter being sampled without taking into account the magnetic field. 

However, other studies have shown that the electromagnetic background can 
have, via quark loop effects, a strong influence also on the distribution of non- 
Abelian fields, thus requiring an unquenched approach. A considerable contribu- 
tion to magnetic catalysis appears to be related to gluon field modifications l59l : 
the magnetic field clearly affects also gluonic observables, like the plaquette 11631 or 
those typically related to confinement/deconfinement, like the Polyakov loop lf57l 
[62 [63). That claims for a more systematic investigation about the effects of electro- 
magnetic backgrounds on the gluonic sector of QCD and about the possible inter- 
play between magnetic and chromomagnetic backgrounds, which could have many 
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implications both at a theoretical and at a phenomenological level (e.g. for cosmol- 
ogy or heavy ion collisions). 

The review is organized as follows. In Sec.|2]we discuss the formulation of QCD 
in external fields, considering both the case of an electromagnetic and of a chromo- 
magnetic field. In Sec. [3]we review studies regarding vacuum properties in external 
fields, in particular concerning magnetic catalysis. Sec.|4]is devoted to a discussion 
about the influence of magnetic or chromomagnetic background fields on deconfine- 
ment and chiral symmetry restoration. In Sec. [5] we extend the discussion about the 
influence that electromagnetic fields may have on the gluonic sector, with a particu- 
lar emphasys on symmetries and on suggestions for future studies on the subject. 



2 Background Fields on the Lattice 

In this Section we will briefly review the methods which are commonly adopted for 
a numerical study of QCD in presence of background fields. There is a difference 
between the case of electromagnetic and chromomagnetic backgrounds: in the for- 
mer, the external field simply adds new U (1) degrees of freedom which contribute to 
quark propagation but are decoupled from the SU (3) degrees of freedom and from 
their updating during the Monte-Carlo simulation; in the latter, the background field 
is made up of the same degrees of freedom which are dynamically updated during 
the simulation, i.e. the SU (3) gauge variables, and one has to follow a Schrodinger 
functional approach. 



2.1 Electromagnetic Fields 

A quark field q propagating in the background of a non-abelian SU (3) gauge field 
plus an electromagnetic field is described, in the Euclidean space-time, by a La- 
grangean \\f(Jp+m)\if, where the covariant derivative 

D ll =d il +igA a ll T a +iqa ll (1) 

contains contributions from the non-Abelian gauge field = A a ^T" and from the 
Abelian gauge field an. Here T" are the St/ (3) generators, g is the SI/ (3) gauge 
coupling and q is the quark electric charge. 

Going to a discrete lattice formulation, SU (3) gauge invariance is naturally pre- 
served Wf\ by requiring that quarks pick up an appropriate non-Abelian phase 
when hopping from one site of the lattice to the other, so that the gauge field 
A^ is substituted by the elementary parallel transports U^(n), which reduce to 
U^(n) ~ 1 + igaAn(n) as the lattice spacing a vanishes in the continuum limit. 

A gauge invariant quantity involving an antiquark at site n and a quark at site 
n + fl is therefore y/(n)U^(n)y/(n + fl) (fl is a unit vector on the lattice). The 
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same approach can be taken for the Abelian field, so that quarks going from site 
n to site n + jU pick up also the Abelian phase Uu(n) = exp(ig /™ dx^a^) ~ 
exp(ieq r a^(n)). A possible discretization of the covariant derivative is then 

D^y-t ^- a (t/ M (n)« M (n)V / (n + A) -^(«-A)4( n_ A)V / (n-A)) (2) 

Of course, in presence of different quark species, carrying different electric charges 
q, the U(\) part of the covariant derivative, m^(h), will change from quark to quark. 
Therefore, as usual, the discretized version of the fermion action is a bilinear form 
in the fermion fields, \fr(i)Mjj v(j), however the elements of M now belong to t/(3) 
instead of 517(3). 

2.1.1 Magnetic Fields 

To take an explicit example, which is directly linked to most of the results that we 
discuss in the following, let us consider the finite temperature partition function 
of two degenerate quarks, e.g. the u and d quarks, in presence of a constant and 
uniform magnetic background field. We shall consider a rooted staggered fermion 
discretized version of the theory, in which each quark is described by the fourth root 
of the determinant of a staggered fermion matrix: 

Z(T,B) = J ^Ue~ SG detM^[B,q u ]detM^[B,q d ] (3) 

where in the standard formulation 

1 4 

Mjj[B,q} =am8jj + - £ T] v {i){u v {B,q)(i) U v (i)8 u ^ 
Z v=l 

- u* v (B,q)(i-v) ul(i-v)8 L j+v) ■ (4) 

is the functional integration over the non-Abelian gauge link variables, So is 
the discretized pure gauge action and u(B,q)i, v are the Abelian links correspond- 
ing to the background field; ; and j refer to lattice sites and TJ !)V are the staggered 
phases. We shall consider two different charges for the two flavors, q„ = 2|e|/3 and 
q c i = — |e|/3. Periodic or antiperiodic boundary conditions along the Euclidean time 
direction must be taken, respectively for gauge or fermion fields, in order to define 
a thermal theory: the temperature T corresponds to the inverse temporal extension. 

A special discussion must be devoted to the issue of spatial boundary condi- 
tions. It is well known that, for a given lattice size (which is usually constrained by 
the computational power available) periodic boundary conditions in space (for all 
fields) are the one which best approximate the infinite volume limit, and are there- 
fore the standard choice in lattice simulations. On the other hand, they place some 
constraints on the possible magnetic fields [68 69 701 I7D . 
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Let us consider a magnetic field directed along the z axis of a three-dimensional 
torus, B = Bz, and let l x and L be the torus extensions in the orthogonal directions. 
The circulation of along any closed path in the x — y plane, enclosing a region of 
area A, is given by Stokes theorem 

a^dxpt = AB (5) 

On the other hand it is ambiguous, on a closed surface like a torus, to state which 
is the region enclosed by a given path: we can choose the complementary region of 
area l x l y — A and equally state 

j> a^dxn = (A - l x l y )B . (6) 

The ambiguity can be resolved either by admitting discontinuities for somewhere 
on the torus, or by covering the torus with various patches corresponding to different 
gauge choices. However it is essential to guarantee that the ambiguity be not visible 
by any charged particle moving on the torus, i.e. the phase factor picked up by a 
particle of charge q moving along the given path must be a well defined quantity. 
Such requirement leads to magnetic field quantization, indeed it is satisfied only if 

exp (\qBA) = exp (iqB(A - l x l y )) (7) 

i.e. if 

9 «=^ (8) 

Lxly 

where b is an integer. It is interesting to notice that this is exactly the same argument 
leading to Dirac quantization of the magnetic monopole charge (in that case the 
closed surface is any sphere centered around the monopole). The quantization rule 
depends on the electric charges of the particles moving on the torus Q in our case it 
is set by the d quark, which brings the smallest charge unit qj = — |e|/3, hence 

I _ 6nb _ 6nb ^ 

Ixly Cp'LxLy 

where L x and L y are the dimensionless lattice extensions in the x, y directions. 

A possible choice for the gauge links, corresponding to the continuum gauge 
field 

a y =Bx; fl M =0 for [i=x,z,t. (10) 

is the following: 



1 Of course we assume all particles living on the torus to carry integer multiples of some ele- 
mentary electric charge, otherwise a consistent quantization of the magnetic field would not be 
possible. 
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u y (B,q)(n)=e ia2 i Bn *; u^(B,q)(n) = I for [i=x,z,t . (11) 

The smoothness of the background field across the boundaries and the gauge invari- 
ance of the fermion action are guaranteed if appropriate boundary conditions are 
taken for the fermion fields along the x direction 1711 . that corresponds to modify- 
ing the 1/(1) gauge links in the x direction as follows: 

u x{B,l)(n)\n x =L x = e~ u,lqLxBny . (12) 

A different condition on the possible magnetic field values explorable on the 
lattice is placed by the ultraviolet (UV) cutoff, i.e. by discretization itself. All the 
information about the presence of the magnetic field is contained in the phase factors 
picked up by particles moving on the lattice. There is a minimum such path on a 
cubic lattice, the plaquette; in particular, a particle moving around ax — y plaquette 
takes a phase factor 

(2Kb \ 
i I (13) 
L X Ly J 

and all other possible phase factors, corresponding to different paths, are positive 
integer powers of that. It is evident that the phase factor in Eq. ( fl"3l cannot distin- 
guish magnetic fields such that qa 2 B differs by multiples of 2k: one would need 
smaller paths to do that, which are unavailable because of the UV cutoff. Hence it 
is possible to define a sort of first Brillouin zone for the magnetic field 

K K 

-—<qB< — (14) 
a L a L 

and we expect that physical quantities must be periodic in qB with a period 2n/a 2 ; 
such periodicity is unphysical and one should always be cautious when exploring 
magnetic field values which get close to the limiting values in Eq. (Tufl i. 

Another issue which is relevant to the lattice implementation of background 
fields regards translational invariance. The magnetic field B breaks the translational 
invariance of the continuum torus in the x — y directions explicitly llTD . That is 
clearly seen by looking at the 1/(1) Wilson lines (holonomies) W x and W y , i.e. the 
t/(l) parallel transports along straight paths in the x or y direction and closed around 
the torus by periodic boundary conditions, which are equal to 

W x = exp(-iqBl x y) ; W y = exp(iqBl y x) (15) 

as can be verified explicitly by means of the elementary parallel transports given 
in Eq. (fTTT i and (fT2l . The Wilson lines therefore leave only a residual simmetry 
corresponding to discrete translations by multiples of 

2k l x . 2k l y 
ax= Wy = b'' U y = qW x = b'' (16) 
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note that Wilson lines are gauge invariant quantities, hence the breaking of transla- 
tional invariance is not a matter of gauge choice. 

However, on a discrete cubic lattice, translational invariance is already broken to 
a discrete residual group, corresponding to multiples of the lattice spacing a, hence 
for a lattice theory in presence of a constant magnetic field, translational invariance 
is further reduced to discrete steps which are multiples of both a x i y and a. Since 
a/a x / y = bjL x j y , lattice translational symmetry is preserved for b multiple of L x / y , 
but is strongly reduced or completely lost for different values, despite the fact that 
the magnetic background field is uniform. 

We have discussed so far about one possible lattice implementation of the 
fermionic action in presence of a magnetic field, Eqs. and ©, which is based on 
the simplest symmetric discretization of the covariant derivative, which takes quarks 
from nearest neighbours sites, rotating them by a non-Abelian and an Abelian 
phase, both corresponding to the simplest straight path, i.e. the elementary paral- 
lel transport. Such kind of discretization has been adopted in many lattice studies, 



It is well known, however, that different, improved discretizations can lead to 
an improved convergence towards the continuum limit. Improvement can proceed 
in different ways: one can consider less simple discretizations of the derivative, 
involving fermion fields which are not nearest neighbours; one can also consider and 
average over paths different from the simplest straight one, so as to smear out UV 
fluctuations: this is idea adopted, for instance, in |62 65 66], where stout smearing 
has been used. If one considers lattice derivatives involving non-nearest neighbours 
lattice sites, one must include appropriate composite U(\) parallel transports, so as 
to preserve the U(l ) gauge invariance of the fermion action. Which paths, however, 
and how the chosen paths are related to the paths considered for the non-Abelian 
phases, is a question which leaves much freedom. 

Let us consider, for simplicity, a gauge invariant bilinear term involving nearest 
neighbour quarks, \j/(n)U^(n)\j/(n + £l), which reduces to ifr(x)(l + aD lx )\j/(x) + 
0{a 2 ) in the naive continuum limit. We can choose a more general, improved such 
term, involving different paths, i.e. we can modify 



where C(«,n + fl) are a set of paths connecting n to n + jl (e.g. for smearing one 
includes staples), Uq are the parallel transports taken along those paths and the co- 
efficients (Xc are chosen in order to keep the naive continuum limit unchanged, up 
to 0(a 2 ). Let us now consider the inclusion of a U(l) external field: a possible 
prescription is that, as quarks explore all considered paths C, taking the associated 
non-Abelian phases, they also take the corresponding U(l) phases along the same 
paths, i.e. 



see e.g. Il57ll59ll63l 




(17) 
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y(n)Un(n)un(n)y(n + p.) -> f(n)l £ OcU c uc\y{n + p) . (18) 

\C(n,«+A) / 

However, different prescriptions can be taken and, in general, the sum over 1/(1) 
phases can be decoupled from that on SU (3) phases, without affecting gauge invari- 
ance. We can consider for example 



where C' is a different set of paths. This is the choice made, for instance, in J62j 
|65]|66), where stout smearing is applied to SI/ (3) links, while the 1/(1) phases are 
left unchanged (i.e. C' runs over the elementary straight path only); that seems a 
reasonable choice, since the 1/(1) field is already a classical smooth field. 

It can be easily shown that for a uniform background and for free fermions, 
i.e. in absence of non-Abelian fields, all possible different choices of paths for 
the 1/(1) term are equivalent, apart from a global phase stemming from interfer- 
ence (Aharonov-Bohm like) effects among the different paths. In presence of non- 
Abelian fields, however, this is in general not true and it would be interesting to 
explore the systematics connected to different choices. 

A final comment concerns the issue of possible renormalizations related to the in- 
troduction of the external magnetic field. The quantity entering the covariant deriva- 
tive and quark dynamics is the combination qa^, which does not renormalize. The 
situation is clearer in the lattice formulation, where the objects carrying informa- 
tion about the external field are the parallel transports of the 1/(1) gauge field along 
closed loops (think of the loop expansion of the fermionic determinant). Such phase 
factors are gauge inviariant quantities and cannot renormalize. The only remaining 
possibility is that the lattice spacing a itself depends on the magnetic field B, be- 
cause of external field effects at the scale of the UV cutoff, so that the physical area 
enclosed by a given loop changes, leading to an effective external field renormal- 
ization. Such possibility, however, has been excluded by a detailed study presented 
in ll62l . showing that a does not depend on B within errors. 

2.1.2 Electric Fields 

After the inclusion of a background magnetic field, the Ip operator is still anti- 
Hermitean and it still anticommutes with 75, hence its spectrum is purely imaginary 
with non-zero eigenvalues coming in conjugate pairs, so that det(lp + m) > and 
Monte-Carlo simulations are feasible, i.e. the path integral measure can be inter- 
preted as a probability distribution function over gauge configurations. 

The situation is different in presence of background electric fields ll72l l42l |4D . 
Let us consider a constant electric field E = Ez, there are various possible choices 
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for a corresponding vector potential in Minkowski space, like (A, = —Ez, A, = 0) 
or (A- = — Et, A, = A y = A, = 0). In every case, after continuation to Euclidean 
space time, t — >• — ix and A, — > iAo, the vector potential becomes purely imaginary, 
thus destroying the anti-Hermitean properties of the Dirac matrix. 

That means that the fermion determinant is complex and numerical simulations 
are not feasible. This "sign problem" has a very strict relation with the usual sign 
problem which is encountered in the study of QCD at finite baryon chemical poten- 
tial, indeed also a baryon chemical potential can be viewed as a constant background 
potential Ao ^ 0. As for finite density QCD, also in this case a possible way out is to 
consider a purely imaginary electric field: the lattice formulation is then completely 
analogous to that of a magnetic field. For example, a purely imaginary electric field 
along the z direction can be formulated on the lattice following Eqs. (fTTT i and (TT2l) . 
by just replacing replacing (x,y) — > (z,t) and B — > Iir^is). A similar approach, fol- 
lowed by analytic continuation to real electric fields, is that usually adopted for the 
study of hadron electric polarizabilities. 



2.2 Chromomagnetic Background Fields 

As we have discussed above, the implementation of an electromagnetic background 
field amounts to adding extra 1/(1) degrees of freedom to the Dirac matrix: the 
added 1/(1) field is a classic, static field which is decoupled from the SI/ (3) gauge 
field appearing in the functional integration. The situation is different for a color 
background field, since in this case one would like to consider quantum fluctuations 
of the non-Abelian field around a given backgroud, e.g. a static ad uniform chromo- 
magnetic field, however the functional integration over gauge variables can destroy 
information about the background field completely. 

A common way to approach the problem is based on the Schrodinger functional 
approach [43 44 1: one considers functional integration over Euclidean space-time 
with fixed temporal boundaries, %\ and x%, and over gauge configurations which 
are frozen to particular assigned values, Af ttl (x,Ti) and Af a2 (x, T2), at the initial 
and final times. This is related to the quantum amplitude of passing from the field 
eigenstate |A? ATl ) to the field eigenstate |A?-° 2 ) in the time {%i — Ti). The amplitude is 
dominated, in the classical limit, by the gauge configuration which has the minimal 
action among those with the given boundary conditions; functional integration can 
then be viewed as an integration over quantum fluctuations around this classical 
background field. 

In Refs. |45l [46] j47l 148] [491 such formalism has been implemented on the lattice 
by considering equal initial and final fields, Af* 1 (x, Ti) = Af ,t2 (x, T2) = A- Xt (x). In 
this way one can define a lattice gauge invariant effective action F[A ext ] for the 
external background field A ext : 




(20) 
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where L, is the lattice temporal extension and i2°[A ext ] is the lattice functional inte- 
gral 



where Sg and Sf are the pure gauge and fermion action respectively, M is the 
fermionic matrix. The functional integration is performed over the lattice links, with 
the constraint 



where U£ xX (x) are the elementary parallel transports corresponding to the continuum 
gauge potential A ext (x) = A^ xt (x)A,„/2; fermion fields as well as temporal links are 
left unconstrained. JF[0] is defined as in Eq. ( I2TI 1. but adopting a zero external field. 

In the following we shall consider the case of an Abelian static and uniform chro- 
momagnetic background field: in this case one can choose X a belonging to the Car- 
tan subalgebra of the gauge group, e.g. X a = A3, while the explicit form of A„ xt (x), 
hence of the lattice links, can be chosen in the same way as for a standard magnetic 
field, see Eqs. ( TTOl ). ( fTTT i and ( fT2b . In this case, as well as in other cases in which the 
static background field does not vanish at spatial infinity, one usually imposes also 
that spatial links exiting from sites belonging to the spatial boundaries are fixed, at 
all times, according to Eq. d22l i: that corresponds to the requirement that fluctuations 
over the background field vanish at infinity. 

Notice that a colored background field influences directly, and through the same 
gauge coupling g, the dynamics of both quark and gluon fields. Indeed numerical 
simulations in presence of colored background fields have considered originally just 
the case of pure Yang-Mills theories. 

Differently from the usual formulation of the lattice Schrodinger functional l43l 
l44l . where a cylindrical geometry is adopted, the effective action defined by Eq. ( f20l > 
assumes an hypertoroidal geometry, i.e. the first and the last time slice are identi- 
fied and periodic boundary conditions are assumed in the time direction for gluon 
fields. For finite values of L t , having adopted also the prescription of anti-periodic 
boundary conditions in time direction for quark fields, Eq. (I2TI 1 can be interpreted 
as the thermal partition function [A ext ] in presence of the background field A ext , 
with the temperature given by T = \/(aL t ). The gauge invariant effective action, 
Eq. ( T20t . is then replaced by the free energy functional: 




(21) 



U k (x,x t = 0) = U k (x,x t = L t ) = l/f l (x) , (k = 1,2,3) , 



(22) 




(23) 
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3 Vacuum Properties in Backgroud Fieds: Magnetic Catalysis 

One of the most significant effects that a magnetic background field can induce on 
the QCD vacuum, as well as on other systems characterized by the chiral properties 
of fermion fields, is known as magnetic catalysis. It consists in an enhancement of 
chiral symmetry breaking, or spontaneous mass generation, which can be thought 
as related to the dimensional reduction taking place in the dynamics of particles 
moving in a strong external magnetic field (see l35ll for a recent comprehensive 
review). 

The enhancement of chiral symmetry breaking reveals itself in a dependence of 
the chiral condensate (y/y/) on the magnetic field B; there are several predictions 
for the actual functional form (y7y/)(.B), depending on the adopted model, how- 
ever some general features can be summarized as follows. By charge conjugation 
symmetry, (y/y/) must be an even function of B, hence, if the theory is analytic 
around B = 0, one expects that (iffy/) depends quadratically on B in a regime of 
small enough fields. The analyticity assumption is violated in presence of charged 
massless fermions, indeed chiral perturbation theory [19] predicts a linear behavior 
in B if m K = 0; such linear dependence is recovered in computations at m K ^ 0, in 
the limit eB l20l . Different power behaviors in B can be found in other model 

computations. In the following we shall make use of the relative increment of the 
chiral condensate, defined as 

V 1 (y/y/)(fl = 0) 

Early lattice investigations of magnetic catalysis have been done in the quenched 
approximation, both for SU (2) l53l and for SU (3) l58l gauge theories. In this case, 
the chiral condensate is computed by inverting a Dirac operator which contains the 
contribution from the external field, but on gauge configurations sampled in absence 
of dynamical fermion contributions. The outcome is that r(B) oc B for SU(2) ll53l . 
while r(B) °< B v , with v - 1.6, forS£/(3) l58l . 

On the other hand, the inclusion of contributions from dynamical fermions mov- 
ing in the background field may produce significant effects. The chiral condensate 
is related, via the Banks - Casher relation 1731 . to the density p(A) of eigenvalues 
of the Dirac operator around A = 0: (y/y/) = izp{0). Such density can change, as 
a function of B, both because the Dirac operator definition itself changes, and be- 
cause the distribution of gauge configurations, over which the operator is computed, 
is modified by dynamical fermion contributions. Studies of magnetic catalysis in- 
volving full dynamical simulations have been reported in [59, 65), for QCD with 3 
colors, and in l63l for QCD with 2 colors. 

In Ref. ll59l . where a standard staggered fermion and plaquette gauge dis- 
cretization has been adopted, corresponding to a pseudo-Goldstone pion mass 
niji ~ 200 MeV, an attempt has been made to separate the contributions to mag- 
netic catalysis coming from the modification of the Dirac operator from those com- 
ing from the modified gauge field distribution: such contributions have been named 
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Fig. 1 We show the relative increment r(B) (see Eq. i24) ) as a function of \e\B for N/ = 2 QCD 
with ~ 200 MeV. Data are reported separately for the u and d quark condensates, as well as 
for their average, together with a best fit according to the functional dependence inspired by chiral 
perturbation theory 1201 . Figure taken from 1591 . 



"valence" and "dynamical", respectively. One can show that such separation can 
be done consistently in the limit of small external fields. Indeed, let us define the 
quantities 

(W) V U %(B) = J @Un^U,0]K(M- x [ mi B Au/d ]) (25) 

(xj/ W )u/"d( B ) = j ' ®U&[m,U,B}Tr(M- l [m,Q Auld ]) (26) 

and, from those, the average (over flavor) quantities (y/y/)™' and (y/y/)''" 1 ; & is the 
probability distribution for gauge fields (including quark loop effects) and M is the 
fermion matrix. In the first case, one looks at the spectrum of the fermion matrix 
which includes the magnetic field explicitly, but is defined on non-Abelian config- 
urations sampled with a measure !0U &[m,U ,0] which, even including dynamical 
fermion contributions, is taken at B = (partial quenching). In the second case the 
measure term takes into account the external field, which is however neglected in the 
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definition of the fermion matrix. From (yy) val (B) and (\jf\jf) < yn (B) we can define 
the corresponding quantities r val ' dyn '. 

In the limit of small fields, B acts as a perturbation for both the measure term 
& [m, U, B] and the observable Tr (M~ 1 [m, B, q u / c i\ ) : assuming quadratic corrections 
in B (however this is not an essential assumption), one can write, configuration by 
configuration: 

0>[m,U,B] = 0>[}n,U,O]+CB 2 + O(B 4 ) (27) 

and 

Tr(M _1 [B]) =Tr(M" 1 [0])+C'B 2 + O(B 4 ) (28) 

where C and C' depend in general on the quark mass and on the chosen configura- 
tion. Putting together the two expansions, one obtains ||59l 

r{B) = r ml (B) + r dy " (B) + 0(B 4 ) , (29) 

therefore, at least in the limit of small fields, the separation of magnetic catalysis in 
a valence and a dynamical contribution is a well defined concept. 

In Fig. [T]we report data obtained in |59l for the relative increment of the u and 
d quark condensates and for their average, as a function of \e\B. Magnetic catalysis 
is quite larger for the u quark with respect to the d quark: this is expected on the 
basis of the larger electric charge of the u quark. Regarding their average, one ob- 
serves that a functional dependence inspired by chiral perturbation theory ll20l fits 
well data, when taking into account the unphysical quark mass spectrum considered 
in l59l . A simple quadratic dependence describes well data in the regime of small 
background fields, approaching a linear behavior for larger fields, while saturation 
effects starts to be visible as x/JcjB reaches the scale of the UV cutoff, which is 
quite low for the lattice setup used in l59l . a -1 ~ 0.7 GeV. Regarding the separation 
into valence and dynamical contributions, Fig. [2 taken again from |59l , shows that 
it is indeed well defined in a significant range of magnetic fields, where the two 
contributions are roughly additive, in the sense that their sum gives back the full 
signal, as expected at least in the limit of small fields. Moreover, from Fig. [2] one 
learns that the dynamical contribution is roughly 40% of the total signal, at least 
for the discretization and quark mass spectrum adopted in 1591 : this is an important 
contribution, which is larger than other usual unquenching effects, which are typi- 
cally of the order of 20%, and reflects a significant modification in the distribution 
of gauge fields, induced by the magnetic background field, which should be further 
investigated by future studies. 

A new investigation of magnetic catalysis for the QCD vacuum has been done 
recently ||65l . making use of an improved gauge and (stout) rooted staggered dis- 
cretization of the theory: results for the average increment of the u and d quark 
condensates have been reported for Nf = 2 + 1 flavors (i.e. including also strange 
quark contributions), adopting physical quark masses and after extrapolation to the 
continuum limit. Also in this case one observes, for the relative increment r(B) (see 
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Fig. 2 Relative increment of the average of the u and d quark condensates as a function of the 
magnetic field. r(B),r val (B),r dyn (B) and r v "' (B) + S ly " (B) are reported separately. From |59l . 



Fig. 1 of Ref. l65l ). a quadratic dependence on \e\B for small external fields, fol- 
lowed by an almost linear dependence as \e\B m\; a nice agreement is found with 
predictions from chiral perturbation theory and from PNJL models (741 . at least for 
not too large fields. It is interesting to notice that, in the regime of small fields, the 
results of 15911 and those of 1651 are compatible with each other if they are rescaled 
by a factor m 2 K (which is different in the two studiesjj the prediction from chiral 
perturbation theory is indeed, in the limit of small fields 



The authors of l63l , instead, have presented an investigation of magnetic cataly- 
sis for QCD with two colors and 4 staggered flavors, which are degenerate both in 
mass and electric charge: that permits, contrary to Refs. [59 65], to avoid rooting 
and the possible systematic effects related to it. The lattice discretization is stan- 
dard, as in (59 1, with also a similar range of pseudo-Goldstone pion masses. Results 
are qualitatively similar to those obtained in Il59ll65l . however the authors try also 



2 The author thanks G. Endrodi for giving him access to the continuum extrapolated data of 1 65 1. 
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an extrapolation to the chiral limit, thus verifying the prediction Ifl9l for a linear 
dependence on \e\B in such limit. 

Till now we have discussed about magnetic catalysis in the QCD vacuum, i.e. 
in the low temperature, deconfined phase, where a consistent picture emerges, both 
from model and lattice computations. The situation is less clear as one approaches 
the high temperature, deconfined regime: results from l62l l65l have shown that the 
growth of the chiral condensate with the magnetic field may stop and even turn 
into an inverse magnetic catalysis at high enough temperatures; such effect may 
have a possible interpretation based on the effects of a strong magnetic field on 
gluodynamics (see I36, l35ll75ll and the discussion in Section 0. Inverse magnetic 
catalysis has been predicted by some model computation 1371 . however in a regime 
of low temperature and high baryon density, which is different from that explored 
in l62ll65l ; recently a possible explanation has been proposed according to which 
inverse catalysis derives from dimensional reduction induced by the magnetic field 
on neutral pions ll76ll . On the other hand, the lattice studies reported in 115711631 do 
not show such effect, reporting instead for standard magnetic catalysis at all ex- 
plored temperatures. We will comment on the possible origin of such discrepancies 
in the following Section, where we discuss about the effects of strong external fields 
on the QCD phase diagram. 

Finally, among the studies aimed at clarifying the effects of strong magnetic 
fields on chiral dynamics, we mention those addressing the determination of the 
magnetic susceptibility of the chiral condensate, both in quenched l54l l58l and in 
unquenched QCD ll66ll . which is part of the total contribution to the magnetic sus- 
ceptibility of the QCD vacuum. In particular, recent unquenched results [66 1 show 
that such contribution is of diamagnetic nature. 



4 QCD Phase diagram in External Fields 

Similarly to what happens with other external parameters, like a finite baryon chemi- 
cal potential, the introduction of a background field, either magnetic of chromomag- 
netic, can modify the phase structure of QCD. The interest in such issue is theoret- 
ical, on one side, since a background field can be viewed as yet another parameter 
of the QCD phase diagram, which can help in getting a deeper understanding of 
the dynamics underlying deconfinement and chiral symmetry restoration. There is 
however also great phenomenological interest, since strong background fields may 
be relevant to the cosmological QCD transition, to heavy ion collisions and to the 
physics of some compact astrophysical objects. The main questions regarding the 
QCD transition that one would like to answer can be summarized as follows: 

1) do deconfinement and chiral symmetry restoration remain entangled, or is a 
strong enough magnetic field capable of splitting the two transitions? 

2) does the (pseudo)critical temperature depend on the background field strength, 
and how? 

3) does the nature of the transition depend on the background field strength? 
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Fig. 3 Behavior of the chiral condensate and of the Polyakov loop, as a function of the inverse 
gauge coupling /3 and for different magnetic field quanta b, on a 16 3 x 4 lattice and for a pion mass 
m % ~ 200 MeV. Figure taken from [57]. 



4) does any new, unexpected phase of strongly interacting matter emerge for strong 
enough background fields? 

Many model computations exist, which try to answer those questions 1177117811791 
[80l[in][H2[I2[M[13][83^ in the following 

we will discuss results based on lattice simulations. The focus will be on aspects re- 
garding deconfinement and chiral symmetry restoration, leaving aside other issues, 
like the possible emergence of a superconductive phase for strong enough magnetic 
fields i3|60]|98]|93. 



4.1 Deconfinement Transition in a Strong Magnetic Background 

Investigating the effects of a background magnetic field on thermodynamics and on 
the phase diagram of QCD necessarily requires an approach in which the presence 
of the magnetic field is taken into account at the level of dynamical fermions. 

A first study along these lines has been presented in 1571 , where finite temper- 
ature QCD with two degenerate flavors has been simulated in presence of a con- 
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stant, uniform magnetic background. A standard rooted staggered discretization (see 
Eq. (O) and a plaquette gauge action have been adopted, with a lattice spacing of the 
order of 0.3 fm, and different quark masses, corresponding to a pseudo-Goldstone 
pion mass ranging from 200 MeV to 480 MeV, with magnetic fields going up to 
|e|fi~0.75 GeV 2 . 

Some results from ll57l are reported in Figs. [3] @1 [5] and [6] In particular, Fig. [3] 
shows the behavior of the chiral condensate and of the Polyakov loop as a func- 
tion of the inverse gauge coupling j3, for the lowest pion mass explored {m K ~ 200 
MeV) and for various magnetic fields (expressed in units b of the minimum quantum 
allowed by the periodic boundary conditions). Remember that the physical temper- 
ature, T = 1/L,a(j3), is a monotonic, increasing function of /3, and that, on the 
16 3 x 4 lattice explored in 1571 , the magnetic field is given by \e\B=(3nT 2 /8)b~ 
(0.03 GeV 2 )/? around the transition. 

Three facts are evident from Fig. [3] The chiral condensate increases, as a function 
of B, for all explored temperatures. The inflection points of the chiral condensate 
and of the Polyakov loop, signalling the location of the pseudo-critical temperature, 
move together towards higher temperatures as the magnetic field increases, mean- 
ing that deconfinement and chiral symmetry restoration do not disentangle, at least 
in the explored range of external fields. The drop (rise) of the chiral condensate 
(Polyakov loop) at the transition seems sharper and sharper as the magnetic field 
increases, meaning that the (pseudo)transition is strengthening. 

Such facts are confirmed by Figs.@]and[5] where the susceptibilities of the chiral 
condensate and of the Polyakov loop are reported, for different values of the mag- 
netic field: the peaks move to higher j3 (hence to higher T) and become sharper and 
sharper as B increases, i.e. their height increases and their width decreases. The in- 
creased strength of the transition is also appreciable from the plaquette (pure gauge 
action) distribution at the pseudocritical coupling, which is reported in Fig. [6] and 
which seems to evolve towards a double peak distribution, typical of a first order 
transition, as B increases. Such hints for a change in the nature of the transition, 
however, have not yet been confirmed by simulations on larger lattices. Regarding 
the increase in pseudo-critical temperature, we notice that it is quite modest in mag- 
nitude and of the order of 2% at \e\B ~ 1 GeV 2 . 

After the exploratory study of Ref. 1571 . new studies have appeared in the lit- 
erature, which have added essential information and also opened new interesting 
questions. Let us first consider the investigation reported in l62l . There are three 
essential differences, with respect to l57l . in the lattice discretization and approxi- 
mation of QCD: 

1) still within a rooted staggered fermion formulation, improved gauge (tree level 
Symanzik improved) and fermionic (stout link) discretizations have been imple- 
mented, and different lattice spacings have been explored, in order to get control 
over the continuum limit; 

2) the authors have explored Nf = 2+1 QCD, i.e. they have considered strange 
quark effects; 

3) a physical quark mass spectrum has been adopted. 
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Fig. 4 Disconnected susceptibility of the chiral condensate as a function of the inverse gauge 
coupling [i and for different magnetic field quanta b, on a 16 3 x 4 lattice and for a pion mass 
m, ~ 200 MeV. Figure taken from [57]. 



Concerning results, instead, the most striking difference with respect to ll57l is 
that the pseudocritical temperature decreases, instead of increasing, in presence of 
the magnetic background, in particular it becomes 10-20% lower for \e\B ~ 1 GeV 2 . 
The authors of 1621 have put the decrease of the critical temperature in connection 
with another unexpected phenomenon that they observe, namely inverse magnetic 
catalysis, i.e. the fact that, at high enough temperatures, the chiral condensate starts 
decreasing, instead of increasing, as a function of B. A possibility pointed out in [65 1 
is that inverse catalysis may originate from the gluonic sector, i.e. the distribution 
of gluon fields may change, as an indirect effect of the magnetic field mediated 
by quark loops, so as to destroy the chiral condensate. We point out that the study 
of ll59l has shown instead that, at least at zero temperature, the modified distribution 
of gluon fields contributes to increase magnetic catalysis (dynamical contribution), 
however the situation may indeed be quite different around and above the decon- 
finement temperature. 

Other aspects pointed out in lf57l have instead been confirmed by 11621 : in par- 
ticular the absence of a clear splitting between deconfinement and chiral symmetry 
restoration, induced by the external field, and an increased strength of the transition 
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Fig. 5 Polyakov loop susceptibility as a function of the inverse gauge coupling /3 and for different 
magnetic field quanta b, on a 16 3 x 4 lattice and for a pion mass m x ~ 200 MeV. Figure taken 
from [57). 



at B 7^ 0, even without evidence for a change of its nature. 



A different approach has been followed by the authors of Ref. [63] : they have 
explored 2 color QCD, making use of a discretization very close to that adopted 
in ll57l . i.e. with the same gauge and staggered action and a similar range of pion 
masses. However, to avoid the possible systematic effects connected with taking the 
square or fourth root of the fermionic determinant, they have considered a theory 
with four degenerate flavors, Nf = 4, all carrying the same mass and electric charge. 
Their results are in quite good agreement with those of ll57l : deconfinement and 
chiral symmetry restoration do not split; the transition gets sharper; the pseudo- 
critical temperature increases as a function of B and no inverse catalysis is observed, 
i.e. the chiral condensate is an increasing function of B for all explored temperatures. 
The increase in T c is also larger than what observed in l57l . being of the order of 
10% at \e\B ~ GeV 2 : that can possibly be explained by the fact that the authors of 
Ref. Il63l make use of 4 flavors, instead of 2, moreover all carrying the charge of the 
u quark, hence the influence of the magnetic field on the system can be larger. 

Another interesting aspect, explored in J63), regards the actual fate of chiral sym- 
metry above the transition and in presence of the background field. The authors have 
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Fig. 6 Plaquette distribution at the pseudocritical temperature for m% ~ 200 MeV and for various 
different values of the magnetic field. Figure taken from 1571 . 



explicitly verified, by performing an extrapolation to the chiral limit, that the chiral 
condensate vanishes in the high temperature phase, i.e. that chiral symmetry is exact 
above T c , also in presence of the magnetic backgroud. It is interesting to notice that, 
in the system explored in [63], the magnetic field does not break any of the flavor 
symmetries of the theory, since all quarks carry the same electric charge. However 
also in the standard case, in which quarks carry different charges, there are still di- 
agonal chiral flavor generators which commute with the electric charge matrix and 
are not broken by the introduction of an external electromagnetic background field, 
so that an investigation similar to that performed in 11631 should be done, to enquire 
about the actual realization of the unbroken symmetries above T c . 

To summarize, present lattice investigations about the influence of a strong mag- 
netic field on the QCD phase diagram have given consistent indications about two 
facts: a magnetic field increases the strength of the QCD transition, however it does 
not change its nature, at least for \e\B up to 1 GeV 2 ; in the same range of external 
fields, no significant splitting of deconfinement and chiral symmetry restoration has 
been detected. 

There is a controversial issue, instead, regarding the location of the pseudocrit- 
ical temperature, T c , which increases as a function of \e\B according to the results 
of 15711631 . while it decreases according to the results of 11621 . It is perfectly possi- 



Lattice QCD Simulations in External Background Fields 



21 



5000 r 



4000 ■ 



3000 



CO. 



2000 



1000 ■ 



P 



10 



Fig. 7 Peak of the derivative of the free energy of the background field with respect to the in- 
verse gauge coupling, together with the chiral condensate, for Nf = 2 QCD and for a background 
chromomagnetic field gB = 7l/(l6a 2 ) ~ 0.35 GeV 2 . Figure taken from |49l . 



ble that various improvements regarding the lattice discretization may turn the very 
slowly increasing function T C (B) determined in l57l into a decreasing function. The 
main point is to understand which aspect is more directly related to the change of 
behavior and to the appearance of inverse catalysis. Is it just a problem of lattice dis- 
cretization and approach to the continuum limit, or can the difference be traced back 
to the larger pion mass used in ||57l ? Could instead the introduction of the strange 
quark explain the differences? We believe that a clear answer to those questions can 
be given by considering each single aspect separately, and that the answer will be 
important by itself, since it will help clarifying the origin of unexpected behaviors 
such as inverse magnetic catalysis. 
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Fig. 8 Peak of the derivative of the free energy of the background field with respect to the inverse 
gauge coupling, the Polyakov loop, for Nf = 2 QCD and for a background chromomagnetic field 
gB = 71 /{16a 2 ) ~ 0.35 GeV 2 . Figure taken from (49). 



4.2 Deconfinement Transition in a Chromomagnetic Background 

A chromomagnetic background field influences directly gluodynamics, hence it 
makes sense to investigate its effects also in pure gauge simulations. A standard 
approach to introduce colored background fields has been described in Section |Z2l 
and is based on the formalism of the lattice Schrodinger functional |43l l44l : such 
approach has been adopted in the literature to study the influence of background 
fields both in pure Yang-Mills theories and in full QCD, and for various kinds of 
background fields, going from those corresponding to a uniform magnetic field to 
those produced by magnetic monopoles [45 46, 47, 48, 49, 50, 51). In the following 
we shall consider the case of a constant background field 147] |48l |49l . 

Like a magnetic field, also a chromomagnetic background field leads to a shift 
of both chiral symmetry restoration, signalled by the drop of the chiral condensate, 
and of deconfinement, signalled by the rise of the Polyakov loop. The two transitions 
move together, moreover it is interesting to notice that the transition point coincides 
with the temperature at which the free energy (effective potential) of the background 
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Fig. 9 Dependence of the critical temperature T c on the external chromomagnetic field for the 
pure gauge theory and for Nf = 2 QCD. Physical quantities are expressed in terms of the parameter 
A ~ 6 MeV. 



field shows a sudden change, the background field being screened (not screened) in 
the confined (deconfined) phase. 

As an example, in Figs. [7] and [8] taken from |49j], we show the behavior of the 
chiral condensate and of the Polyakov loop versus the inverse gauge coupling, j3, to- 
gether with the derivative of the free energy with respect to /3, which shows 
a sharp peak at the same point at which chiral symmetry restoration and decon- 
finement take place. Results make reference to simulations of Nf = 2 QCD on a 
32 3 x 8 lattice, with a standard pure gauge and staggered fermion formulation, a 
lattice spacing a ~ 0.15 fm and a bare quark mass am = 0.075; the value of the 
magnetic background field is gB = %/(\6a 2 ) ~ 0.35 GeV 2 . 

Present lattice results indicate that, both for pure gauge and full QCD, the tran- 
sition temperature decreases in presence of a constant background chromomagnetic 
field. This is shown in Fig. [9] [49], where the critical temperature T c (expressed in 
units of the parameter A ~ 6 MeV), is shown as a function of gB, both for the pure 
gauge theory and Nf = 2 QCD. The change in the critical temperature is, in gen- 
eral, much larger than what happens in presence of an electromagnetic background: 
this can be interpreted in terms of the fact that a colored background field directly 
affects gluodynamics. An extrapolation of lattice results would even hint, as shown 
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Fig. 10 Chiral condensate as a function of the inverse gauge coupling /3 for different values of the 
external chromomagnetic field gB = 7tn ext /(16a 2 ) and for Nf = 2 QCD. Figure taken from |49|. 



in Fig. [9] at the presence of a zero temperature deconfining transition, for y/gB of 
the order of 1 GeV. It would be interesting to investigate this possibility further, 
in the future, for the possible cosmological and astrophysical implications it could 
have. One should repeat the study of 11491 in presence of more physical quark masses 
(with the values adopted in (49), m K is of the order of 400-500 MeV) and closer to 
the continuum limit, also in order to reach higher values of the external field. 

As a final comment, we notice that |49l also reports evidence about magnetic 
catalysis induced by the chromomagnetic field. Fig.[l0] in particular, shows the be- 
havior of the chiral condensate for a few values of gB: it is clear that, at least around 
the transition, where data for all external fields are available, the chiral condensate 
grows as a function of gB. It is interesting to notice that in this case T c decreases 
anyway, even in presence of normal catalysis, i.e. inverse catalysis does not seem to 
be a necessary condition for a decreasing critical temperature, at least in presence 
of finite quark masses. 
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5 More on Gauge Field Modifications in External 
Electromagnetic Fields 

The lattice results that we have discussed in the previous Sections show that an 
electromagnetic background field can have a strong influence not only on quark dy- 
namics, but also on gluodynamics, even if the interaction is not direct but mediated 
by quark loop effects. That does not come completely unexpected, at least in the 
strongly interacting, non-perturbative regime. Let us summarize a few facts: 

1) we have shown that large part of magnetic catalysis is due to the modification of 
the gauge field distribution (dynamical contribution), i.e. that the increment of chiral 
symmetry breaking is lower by about 40% if gauge field configurations are sampled 
without taking into account the background field |59l ; 

2) it is known, from model computations, that a strong electromagnetic background 
can modify gluon screening properties and influence confinement dynamics, reduc- 
ing the confinement scale ll36l l75l l35l : such effects on gluodynamics have been 
proposed in |65l as a possible explanation for inverse catalyis; 

3) lattice computations have also explicitly shown that quantities related to quark 
confinement, like the Polyakov loop, are influenced by a background electromag- 
netic field, not only around the transition but also well inside the confined and the 
deconfined phases (see as an example the data reported in Fig.|2|i; 



4) the results of [63] have shown that the magnetic field can also induce an asymme- 
try in the non-Abelian plaquette values, especially around and above the transition, 
with a possible significant effect on the equation of state of the Quark-Gluon Plasma. 

It is surely of great importance to study and understand these effects in a deeper 
and more systematic way. Is it possible, for instance, that electromagnetic back- 
ground fields may induce chromoelectric or chromomagnetic background fields? 
That has been discussed in some recent model studies (see e.g. 1751 ) and could have 
important phenomenological consequences, both for cosmology and astrophysics 
and for heavy ion collisions. In the following we would like to discuss how lat- 
tice simulations could further contribute to the issue, with an accent on symmetry 
aspects. 

Let us consider a uniform and constant magnetic background field. In order to 
understand some of the gluonic field modifications induced by it, it is interesting to 
observe that it breaks explicitly charge conjugation symmetry, and to ask how such 
breaking propagates to the gluon sector. 

Charge conjugation on gluon fields means — > — or, at the level of lattice 
gauge link variables, C/u — > t/^. Symmetry under it means that a gauge configuration 
Ufi{x) has the same path integral probability of its complex conjugate U£(x): 



@HJ@>[U}=giU*&>[U*]. 



(31) 



As a consequence, all gluonic gauge inviarant quantities, i.e. traces over closed par- 
allel transports Wc[I/], including plaquettes and Wilson loops, must be real: 
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{W C [U})* = (W£[U]) = (W C [U*}) = J Z>U0>[U}W C [U*} 

@U*0>[U*]W c [U] = f @U&>[U]W C [U} = (W C [U\) (32) 



where the first equality holds if the path integral measure is real, the second express 
the fact that Wc is a trace over a product of links, while the fourth is simply a 
change of integration variables. Eq. d32l may be violated by possible spontaneous 
symmetry breaking effects, which may affect some particular loops, as it happens 
for the Polyakov loop in the deconfined phase. 

In presence of an electromagnetic background field, such property must be lost, 
since the breaking of charge conjugation symmetry will propagate from the quark 
to the gluonic sector: £P[U] ^ £P[U*] because of the contribution from the fermion 
determinant. Let us see that more explicitly, considering the particular case of the 
trace over a plaquette. 

The most direct way to look for effective quark contributions to the gauge field 
distribution is to consider the loop expansion of the fermionic determinant. If we 
write the fermion matrix as M = nifld + D, where D is the discretization of the 
Dirac operator and m f is the bare quark mass, then the following formal expansion 
holds: 

detM = exp(Trlog(m / Id + D))ocexp^-^y^]-TrD A ^ . (33) 

D k is made up of parallel transports connecting lattice sites; the trace operator im- 
plies that only closed parallel transports will contribute. Moreover, for each con- 
tributing loop, there is an equal contribution from its hermitian conjugate. If D is a 
standard, nearest neighbour discretization, the first non-trivial term comes for k = 4 
and contains a coupling to the plaquette operator, which can be expressed as: 

MMI (TrH MvW+ Trnt vW ) = MM Re TrII, v (x) (34) 

where Aj5(mf) is a coefficient which depends on nif and on the particular fermion 
discretization adopted. The total effect can be viewed as a simple renormalization 
of the inverse bare gauge coupling, j8 — > j3 + Y,f^P(mf), where the sum runs over 
flavors. 

Let us now consider the effect of an external electromagnetic field. That will 
change the definition of D, see e.g. Eq. (0, so that each loop contribution to the de- 
terminant will get a 1/(1) phase from the external field. In particular, the expression 
in Eq. (l34l for plaquettes will be modified into 

MM I UvW M(lvW + e - i frvw TrJT t v(x; 

Nc 2 V 

= ^. ( . OT/) (cos(0 MV (*))ReTrIT M vOe) - sin(^ v (x))ImTrn MV (x)) . (35) 
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The phases 0^ v are in general non-trivial. For the particular case of a uniform mag- 
netic field B in the z direction we have ^ = qfBa 2 , where qy is the quark charge. 
This implies that, apart from the case in which there is exact cancellation among 
quark flavors carrying equal masses and opposite electric charges (but this does not 
happen in real QCD), there will be a non-trivial coupling to the imaginary part of 
some plaquette traces, which will then develop a non-zero expectation value. No- 
tice that this is exactly what happens also for the Polyakov loop in presence of an 
imaginary baryon chemical potential. 

It is interesting to understand the meaning of a non-zero expectation value for 
the imaginary part of the trace of the plaquette, in terms of continuum quantities. 
Remember that, in the formal continuum limit, the plaquette operator is linked to 
the gauge field strenght G^v = G c ^ v T a as follows: 

n MV ^ exp (i G MV fl 2 ) (36) 

Considering the expansion of the exponential on the right hand side, the first non- 
trivial contribution to the imaginary part of its trace comes from the third order, 
which indeed is proportional to 

-iTrG^ v = -iG; v G%G^Tr{T"T h r) = ^ -irf*) GJ V G* V G£ V (37) 

where no sum over , v is understood. A non-zero expectation value for the imag- 
inary part of the plaquette, therefore, corresponds to a non-zero charge odd, three 
gluon condensate, constructed in terms of the symmetric tensor d abc : 

(d^G^G^G^^O. (38) 

Such condensate, which is zero by charge conjugation symmetry in the normal 
QCD vacuum, will develop a non-zero value along some directions, e.g. for (/i, v) = 
(x,y) in case of a magnetic field in the z direction, and can be considered as a non- 
trivial, charge-odd gluon background induced by the external field. Of course the 
effect will be present only for N c > 3 colors, since for SU (2) all traces are real and 
indeed d" bc = 0. 

The loop expansion of the determinant, examined above, is just a way to show 
explicitly that such effects can be present, but cannot be quantitative, since the ex- 
pansion is purely formal (apart from the limit of large quark masses). A systematic 
investigation can and must be performed in the future by means of lattice simula- 
tions. 

Similarly to charge conjugation symmetry, one can investigate electromagnetic 
backgrounds which break explicitly the symmetry under CP, i.e. such that E • B ^ 0, 
and ask how CP violation propagates to the gluon sector, giving rise to an effective 
parameter. Such phenomenon is in some sense complementary to the chiral mag- 
netic effect and is related in general to the effective QED-QCD interactions in the 
pseudoscalar channel H100I UOll 11021 [T03l . It has been studied recently by a first 
exploratory lattice investigation 01041 . where numerical simulations at imaginary 
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electric fields plus analytic continuation have been exploited to determine the sus- 
ceptibility Xcp of the QCD vacuum to CP-odd electromagnetic fields, defined by 
e ff ~ ^c/>e 2 E ■ B, obtaining Xcp ~ 7 GeV~ 4 for QCD with two staggered flavors 
and m n ~ 480 MeV. 

Finally, let us discuss a few other lines along which the effects of electromagnetic 
background fields on the gluon sector could be further investigated. The first consists 
in the determination of the effective action of given gluonic backgrounds, using the 
techniques described in Section lzSl but in combination with non-trivial electromag- 
netic external fields: that would give the possibility of studying if the latter can give 
rise to instabilities in the gluon sector, leading to the generation of non-trivial glu- 
onic backgrounds (as proposed e.g. in 11751 ). The second regards the measurement of 
glueball masses and screening masses in presence of magnetic backgrounds, which 
could give the possibility, for instance, to test the proposal given in ||361 about the 
lowering of the confinement scale in presence of a strong magnetic field. Finally, it 
would be important to perform a careful investigation about the influence of mag- 
netic backgrounds on the QCD equation of state, as proposed in l63l . 



6 Conclusions 

Numerical simulations in presence of background external fields have been consid- 
ered since the early stages of Lattice QCD. The last few years, however, have seen 
the development of considerable activity on the subject, which has been driven by 
the recent interest in theoretical and phenomenological issues regarding the behav- 
ior of strongly interacting matter in magnetic fields. 

In this review, after a general overview about the formulation of lattice QCD in 
presence of magnetic or chromomagnetic background fields, we have discussed the 
present status of lattice studies regarding magnetic catalysis and the QCD phase 
diagram. In the low temperature phase, consistent results have been obtained, by 
different groups, which confirm magnetic catalysis and the predictions coming from 
chiral perturbation theory and some effective models. Good part of the enhancement 
of chiral symmetry breaking seems to be associated with modifications of gluon 
fields induced by the magnetic field via quark loops. A chromagnetic background 
induces magnetic catalysis as well. 

In the high temperature phase, consistent results have been obtained regarding the 
fact that a magnetic or chromomagnetic field shifts both chiral symmetry restoration 
and deconfinement and that the two transitions do not split, at least for fields up to 
0(1) GeV 2 . A magnetic field also leads to a strengthening of the transition, which 
however does not seem to turn into a strong first order, at least for fields up to 
0(1) GeV 2 ; it would be interesting however, to further investigate the issue also in 
combination with other external parameters, like a baryon chemical potential, since 
that could be relevant to the search for a critical endpoint in the QCD phase diagram. 

Regarding the shift in the pseudocritical temperature induced by a background 
magnetic field, existing studies have shown that a modification of the lattice imple- 
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mentation of QCD can change the outcome, going from a slight increase (of the 
order a few % for \e\B of (9(1) GeV 2 ) to a decrease (of the order of 10 % for \e\B 
of (9(1) GeV 2 ) when an improved discretization and a physical quark mass spec- 
trum are used; in the latter case one also observes inverse magnetic catalysis, i.e. 
a decrease of the chiral condensate as a function of B, in the high temperature, de- 
confined phase. A decrease of the pseudocritical temperature is also observed in 
presence of a chromomagnetic background field. A remaining question is which ef- 
fect is more relevant to explain the discrepancy among different lattice studies, i.e. 
whether lattice artifacts, or the unphysical quark spectrum, or both are at its origin. 

One of the aspects that emerges more clearly from present studies is the fact that 
gluon fields are strongly modified by external electromagnetic fields, even though 
indirectly, by means of dynamical quark loops. We believe that this aspect can and 
should be investigated more sistematically by lattice simulations: in Sec. [5] we have 
discussed and suggested a few possible ways to do that by future studies. 
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